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ABSTRACT 

The secular perturbations exerted by an inclined satellite orbiting in a gap in 
a broad planetary ring tends to excite the inclinations of the nearby ring particles, 
and the ring's self-gravity can allow that disturbance to propagate away in the 
form of a spiral bending wave. The amplitude of this spiral bending wave is 
determined, as well as the wavelength, which shrinks as the waves propagate 
outwards due to the effects of the central planet's oblateness. The excitation 
of these bending waves also damps the satellite's inclination I. This secular I 
damping is also compared to the inclination excitation that is due to the satellite's 
many other vertical resonances in the ring, and the condition for inclination 
damping is determined. The secular / damping is likely responsible for confining 
the orbits of Saturn's two known gap-embedded moons. Pan and Daphnis, to 
the ring plane. 

Subject headings: planets: rings 



1. Introduction 

Secular gravitational perturbations can play a significant role in determining the global 
structure and the long-term evolution of a disk-companion system. A well-known example 
is the circumstellar dust disk at j3 Pictoris, whose broad but gentle warp is thought to 
be due to the secular gravitational perturbations exerted by an unseen planetary system 
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( iMouillet et al.lll997l ). Sec ular perturbations f rom an eccentric planet can also make a dust- 
disk appear lopsided, too (IWyatt et al.lll999l ). Secular perturbations are those forces that 
are due to the time-independent part of a companion's gravitational potential, and such 
perturbations are equivalent to th e gravitational forces that arise when a perturber's mass 
spread about its orbital ellipse ( Murray fc Dermottlll999l ). Consequently, the long-term 



IS 



secular evolution of a disk- companion system is conveniently modeled by treating it as a 
system of gravitating rings (jHahrul2003l ). an approach that will also be employed here. 



When the perturber's orbit is eccentric or inclined, its secular perturbations can excite 
the orbital eccentricities or inclinations of the disk particles, as well as cause the orbits of the 
disk particles to precess over time. Large eccentricities e or inclinations / can also be excited 
at a secular resonance, which is a site where the disk matter precesses in sync with one of 
the eigenfrequencies that describe the perturber's precession. However, substantial e's and 
J's can also be excited elsewhere at non-resonant sites in the disk, with greater excitation 
occurring nearer the perturber. In fact, it is these non-resonant secular perturbations of the 
disk that are the focus of this study. 

If the disk has internal forces, such as pressure or self-gravity, then those internal forces 
can also transmit the companion's disturbances across the disk. For instance, an inclined 
planet orbiting in a circu mstellar gas disk can ex cite a global warp that is facilitated by the 
disk's internal pressure (iLubow &: Ogilvid l200ll ). And if an eccentric companion inhabits 
a gap in the gas disk, its secular perturbations can launch a densi ty wave at the ga.p edg e 



having such a long wavelength that a global standing wave emerges (iGoldreich &: Sari 



20031 ). 



But if the disk is instead gravity-dominated, then the compani on can launch spiral densit y 
or spiral bending waves at its secular resonances in the disk (IWard fc HahnI Il998l . l2003l ). 
These phenomena are also relevant to studies of extra-solar planets, since any dissipation in 
the disk facilitates a transfer of angular momentum between the disk and the companion in 
a manner that tends to drive its eccentricity or inclination to zero. 

The following will examine the secular evolution of a related system: of a small satellite 
that inhabits a narrow gap in a broad planetary ring, both of which are orbiting an oblate 
central planet. It will be shown below that an inclined satellite can launch a spiral bending 
wave that propagates outwards and away from the gap's outer edge. The amphtude and 
wavelength of this spiral bending wave is assessed below, as well as the rate at which this 
wave-action damps the satellite's inclination. This secular inclination-damping mechanism is 
then compared to the inclinat ion excitation that is due to the satellite's many other vertical 
resonances in the ring (e.g., iBorderies et al.l Il984j ). We then quantify when this secular 
damping dominates over the resonant excitation, and show that this secular interaction is 
likely responsible for confining Saturn's two known gap-embedded moons. Pan and Daphnis, 
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to the ring plane. 



2. Equations of motion 



Begin by considering a planetary ring that is perturbed by a single satellite, with both 
orbiting an oblate planet. To assess the disturbance that the satellite might launch in this 
ring, the Lagrange planetary equations will be used; they give the rates at which a ring 
particle's orbital inclination / and longitude of ascending node Q varies with time t: 



1 OR 



na^I dVt 



and VL 



1 dR 
na^I dl ' 



where R is the disturbing function for a ring particle having a semimajor axis a and mean 
motion n ^ y^/GM/a^, where G is the gravitation constant and M is the mass of the central 
planet ( iMurray fc Dermottlll999l ). and all inclinations are small, / ^ 1. The total disturbing 
function for a ring particle is R = Rdisk + Rsat + Robh where the three terms account for the 
gravitational perturbations that are due to ring's gravity (which we treat here as a broad 
disk), the satellite's perturbations, and that due to the planet's oblate figure. The particle's 
equations of motion is thus the sum of three parts: 



disk 



and = 



5 at 



disk 



+ n 



sat 



obi 



(2) 



noting that oblateness does not alter inclinations. And because we are only dealing with the 
syste m's secular perturbations, the semimajor axes of all bodies are constant (IBrouwer fc Clemence 
196lh . 



The amplitude of a spiral bending wave that is in a steady-state does not vary with 
time, so the disk inclinations obey /(a) = throughout the disk. A persistent spiral pattern 
must also rotate with a constant angular velocity cj, so 



disk 

and uj 



sat 



disk 



+ Q 



+ Q 



sat 



constant. 



obi 



(3) 
(4) 



The following subsection will use the first equation to solve for the wave amplitude /(a) 
throughout the disk. The next subsection will then use the other equation to solve for the 
bending waves' dispersion relation uj{k), which in turn provides the wavenumber k of the 
spiral bending wave, and the wave's radial velocity. 
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2.1. wave amplitude 



Begin by examining how the planetary ring perturbs itself. The ring is to be regarded 
as a broad disk that is composed of many narrow, concentric annuli. Each annulus has 
mass Sm{a), inclination /(a), and longitude of ascending node Q{a), all of which are to be 
regarded as functions of the rings' semimajor axes a. For the moment we will assume that 
all rings are cir cular, notin g that we will deal with the system's eccentricity evolution in a 
foUowup study (lHahnll2007l ). Suppose that the annulus at a is perturbed by another annulus 
of mass 6m' and radius a'; the disturbing function for the perturbed annulus is 



SR 



GSm' 

-abr, Lia) 



4a 



^3/2 V 



II' cos{n-n') 



(5) 



where a, /, Q are the orbit elem ents of the perturbed annulus, and the primed quantities 
refer to the perturbing annulus (IHahnl 120031 ) . The softened Laplace coefficient appearing in 
the above is 

cos{(p)d(p 

I TT^r- — — 

TX 



(6) 



,0 [(l + a2)(l + [)2) _2acos(^]3/2' 

and it is a function of the semimajor axis ratio a = a' /a, where 1) = /i/a <^ 1 is the disk's 
vertical scale height h in units of semimajor axis a. Note that when the disk is infinitesimally 
thin, t] = 0, and the disturbing function 6R is equivalent to that due to a point-mass 6m' 
(e.g., iBrouwer fc Clemencdll96ll ). 



It will be convenient to replace the ring mass 6m' with 27Ta'a'da', where a' = a{a') is 
the mass surface density of the annulus of radius a' and radial width da'. We will also write 
its semimajor axis as a' = a(l + x'), where x' = (a' — a)/a = a — 1 is the fractional distance 
between the perturbing ring a' and the perturbed ring a. For the moment we will consider 
a one-sided disk — one that orbits wholly exterior to the satellite, where A is the fractional 
distance between the satellite's orbit and the disk's inner edge; the geometry is sketched in 
Fig. [TJ The disturbing function for ring a due to perturbations from ring a' can now be 
written as 



6R = -l^;^(na 



.X 



-i^-ii' cos{n-n') 



dx' 



(7) 



where fi'^ = 'Ka'a''^/M is the ring's so-called normalized disk mass, dx' = da' /a is the 
perturbing ring's fractional width, and ^g/gl^^O is shorthand for Eqn. dSj) evaluated at a = 
1 + x'. Then according to Eqn. ([1]), ring a' will alter the inclination of ring a at the rate 



61 



1 di6R) 
na^I dVt 



^'j^nby^{x')I' sin{Q — Q')dx'. 
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2.1.1. ring-disk evolution 

The total rate at which the entire disk alters the inclination of ring a is the above with 
x' integrated across the disk, from —x to +cxo (see Fig. [1]), so 



disk 



SI = -n I ^i'^{x'p^j^{x')I'{x') sin(fi - n')dx' . (9) 



disk 



As one might expect, this integral is dominated by the contributions from nearby annuli 
that lie a small distance x' away. In the \x'\ -C 1 limit, the softened Laplace coefficient is 

*^'/.(-') - ^^51^ (10) 



(Hahn 2003 ). Because of the steep dependence of fei/l on x', we can replace the inclination 



^3/2 

I'{x') and disk mass fi'^ix') with their values evaluated at the perturbed ring at x' = 0, so 
r c:^ I and yu'^ ^ = Tiaa^/M, and also pull them out of the integral so that 

disk vr J^^ x'2 + 2P|2 

Most of the contributions to this integral will be due to nearby annuli that lie a wavelength 
A ~ 27r/|fc| away, where k is the wavenumber of the spiral bending wave. 



IS 



A spiral wave has a wavenumber k{a) = —dVt/da (Eqn. IA2p . so the — fi' in the above 

Q{a) - fi'(a') = - / k{r)dr. (12) 



In general, the wavenumber k{a) will vary with semimajor axis a. However, considerable 
progress can be made if we assume that k is constant over a wavelength, so i7 — fi' ~ 
—k{a — a') = kax'. Then Eqn. (fTT!) becomes 



disk TT 



—AH{\k\ax)fidIkan (13) 



after replacing the x' integration variable with y = \k\ax', and noting that Eqn. (ITT!) is 
odd in y. In the above, the function A{z) is a dimensionless measure of the warped disk's 
perturbation of itself: 



oo 



where z = \k\ax is the distance from the ring-edge in units of 27i wavelengths, and the di- 
mensionless wavenumber H = \^[)\k\a is roughly the disk's vertical thickness in wavelength 
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units. The function Ah{z) is also plotted in Fig.[2l We will be interested in a disk whose ver- 
tical thickness is small compared to t he wavelength, so H <^ 1 , and Ah{z) ^ sm{z)/z — Ci{z), 
where Ci{z) is the cosine integral of lAbramowitz &: StegunI (119721 ) . Also keep in mind that 
these results assumed that the wavenumber k varies little over a single wavelength; subsection 
12.2.21 will note when this approximation breaks down. 



2.1.2. ring-satellite evolution 



The ring at semimajor axis a is also being perturbed by the satellite, and that ring's 
disturbing function Rg due to the satellite is Eqn. with 6171' replaced by the satellite's 
mass m,: 



1 



1(1) 



a) 



-P - lis cos(l] - Qs 



(15) 



where fig = mg/M is the satellite's mass in units of the central planet's mass, and a = 
as/a = (1 + A + x)^^ ~ 1 — (A + a;). The satellite's perturbation thus causes the ring's 
inclination to vary at the rate 



sat = = rsnabflia)lssini^ - fi.). 



For a ring in the wave-excitation zone, i.e., near the satellite 

2 

^3/2 



(16) 



(17) 



7r(x + A)2' 

since the ring lies a fractional distance a; + A away from the satellite (see Fig. [1]), with both 
presumably well separated such that A ^ f). We also write the above longitude difference as 
Q — Qs — —kax + 0o, where the angle 0o allows for the possibility that the annulus nearest 
the satellite at x = may have a longitude of ascending node that differs from the satellite's 
node VLs by angle (po- Thus 

- o T'^^A^-^ sm{-kax + 0o) (18) 
sat 27r(x + A)^ 

is the rate at which the satellite alters a ring's inclination. 



2.1.3. wave amplitude 

When the wave is in steady-state, the two /-excitation rates, Eqns. ffT^ and ffT5]l . are 
balanced, which yields the amplitude of the bending wave: 

I{z) \k\a Us sin(z — Sfe0o) 



2 iid{z+\k\a/\YAH{z) 



(19) 
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where Sk = sgn(/c) and z = \k\ax is the downstream distance in units of 27r wavelengths. 
Far downstream, wh ere z ^ 1, we expect I(z) — » constant. For z ^ 1, Ci{z) ^ sm{z)/z — 
cos{z)/z^ + 0{z~^) (lAbramowitz fc StegunI Il972l ). so Ah{z) ~ cos{z)/z^ downstream; see 
Fig. [2J So if I{z) is to be a finite constant, then the longitude offset must be (^o = ±7r/2, 
and I/Is — — (|fc|o./2)(/is//id)sA:Sin0o- Of course, these inclinations must also be positive, so 
-Sk, and the bending wave amphtude becomes 



sm 



±1 



I 

z 



|fco|Q/is 



(20) 



where |A;o| is the initial wavenumber at x = 0, where the wave is excited at the disk's inner 
edge. To make further use of this result, we still need the initial wavenumber k^, which we 
get from the waves' dispersion relation. 



2.2. dispersion relation 



The waves' dispersion relation is obtained from Eqn. (jlj), with each term in that equation 

is the rate at which the disk drives its own precession. 



assessed below. The first term, Q 



disk 



Again we calculate that rate by treating the disk as numerous narrow annuli. The rate that 
annulus a precesses due to the secular perturbations from the annulus at a' is 



6n 



1 d{6R) 
na?I dl 



(1) 

'3/2 



X 



I'(x') 



dx' 



(21) 



disk 



where 5R is Eqn. ([7]). The total precession rate due to the disk's self-gravity is (l 

J^.^j^6Cl, where the integration proceeds across the entire disk. Again, the integrand is a 
steep function of x', due to the softened Laplace coefficient, Eqn. (ITUI) . which allows us to 
replace the quantities I'{x') and ^i'^{x') with the constants / and ^a- The disk's precession 
rate due to its self-gravity then becomes 

sm^{\k\ax' /2) 



disk 



:fJ'dn 



x'2 + 21^2 



-dx' 



where 



-BH{\k\ax)\k\a^dn 
sm\y/2) 



dy. 



(22) 
(23) 

(24) 



The function Bh{z) is a dimensionless measure of the rate at which the disk drives its own 
precession. When the disk is much thinner than the wavelength, H = y/2[)\k\a -C 1, and 

cos z — 1 



Bu{z) + -Si(^) 
Z n 



TCZ 



(25) 
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where Si{z) is the sine integral of lAbramowitz fc StegunI (119721 ) . Far downstream, where 
z oo, the Bh integral evaluates to 



rjoo 



lim Bh{z) = -(1 

2— >00 Si 



(26) 



Note that B'^ is maximal when the dimensionless wavenumber is small, i.e., H <^ 1, for 



which B 



H 



1. But if the disk is thick, H ^ 1 and B 



H 



H ^, which indicates that the 
disk's ability to sustain a bending wave is weakened when the disk is too thick. 

Figure [2] also shows a numerical evaluation of Bh{z) for a thin disk having H = 0.01. 
This Figure shows that Bh{z) takes values of 1/2 < Bh{z) < 1 for ^ > 0, with ^^^(0) = 1/2 
at the disk's inner edge, and that Bh{z) 1 downstream where z ^ 1, provided the disk 
is thin. 

The satellite is also precessing the ring material orbiting nearest it; that precession 
occurs at the rate 



sat 



1 dRs 
na^I dl 



1 — y- COs(i7 — 



_27r(2; + A)^ 



+ 



[Id sin ( A; 


ax) 


7r 


k 


a{x + A)2 



n 



(27) 
(28) 



where Ig/I is replaced by the downstream wave amplitude, Eqn. (120|) . The first term is the 
familiar differential precession that would occur if the disk were massless. The second term, 
which is proportional to the disk mass fid, is the additional precession that is due to the 
torque that the satellite exerts upon the disk's spiral pattern. 

The central planet's oblateness is also driving precession; the disturbing function for 
that perturbation is 



obi 



4 



Rr 



[an] 



(29) 



wher e J2 is the planet's second zonal harmonic, and Rp is the planet's radius (iMurray fc Dermott 
19991 ). Precession due to oblateness is then 



1 dR, 



■obi 



obi 



na^I dl 

l--(x 
2^ 



3 [Rp 



n 



A) 



obi 



(30) 
(31) 



obi 



where Eqn. fl3Tl) is a Taylor expansion of Eqn. fl30|) in the small quantity x + A, and Clg 

— {3J2/2){Rp/asYns is the rate at which the satellite's orbit precesses due to oblateness, 
where is the satellite's mean motion. 
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Summing Eqns. fl23l) . fl28|l . and fl3Tl) provides the dispersion relation for the spiral bend- 
ing waves: 



l-^(a; + A) 



where 

Sin 2 , , 

^W = ^''W + 7(ITS^- '''' 

All terms in Eqn. fl52]) are negative, so the disk precesses in a retrograde sense. Note that 
if a spiral bending wave is to persist over time, then all parts of the disk must precess in 
concert. The dispersion relation, Eqn. (132|) . thus tells us how the wavenumber \k{x)\ must 
adjust throughout the disk in order for the spiral bending wave to precess coherently. 

The first term in the dispersion relation is due to the disk's self-gravity. That term is 
proportional to D{z), and it has two parts: self-precession that is driven by the bending 
wave itself (the term D), and the additional precession that is driven by the satellite's 
torque on the spiral wave pattern [the latter term in Eqn. fl55]) ]. The function D{z) is plotted 
in Fig. El which shows that 1/2 < D{z) < I. 

The second term in Eqn. (l32l) is the rate at which the satellite drives differential preces- 
sion in the disk; this effect is most prominent nearest the satellite. The third term is the rate 
at which the oblate central planet drives differential precession, and this occurs all through- 
out the disk. Differential precession can inhibit wave-action by shredding the spiral pattern. 
But inspection of the dispersion relation suggests that bending waves can propagate, despite 
differential precession due to the satellite, when the satellite's mass is sufficiently small, i.e., 
when /is ^ yUrf|A;|aA^. The dispersion relation also tells us that the wavenumber \k\ must 
also increase with radial distance x in order to compensate for the additional differential 
precession that is due to the oblate central planet. 



2.2.1. group velocity 



The waves' group velocity is (jToomrelll969l : iBinney &: Tremaind 119871 ) 



dk 



(34) 



upon setting D[z) 1 downs tream; this is the rate at which the spiral bending wave 
propagates radially (lHahnll2003l ). Since the satellite is launching outward-propagating waves 
from the disk's inner edge, the group velocity must be positive, which implies that Sk = 
sgn{k) = —1. Spiral waves having k < are called leading waves. Note also that sin^o = 
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—Sk = +1, so (po = 7r/2, which means that the longitude of ascending node at the disk's 
inner edge leads the satellite's node by 90°. 



2.2.2. wavenumber k 



The wavenumber k can be obtained by calculating the satellite's precession rate Clg. 
When the system is in steady-state, both the satellite and the spiral wave precess at the 
same rate, Qg = ^^{\k\), which provides another equation for the wavenumber k. 



The satellite's node fig is being precessed by the disk and by the central planet, so 



disk 



+ ^s 



obi 



(35) 



where 



disk 



f^._^f^ 6Qs is the satellite's precession rate due to the entire disk, and where 



cos(i7s — Q') 



dx' 



(36) 



is the satellite's precession rate due to a disk annulus of radius a' and mass 6m'. This 
can be obtained from Eqn. (|2Tll with n, a, /, Q replaced by Ug, ag, Is, ^s and the separation 
x' —>■ x' + A. The satellite's precession rate due to the entire disk is 



disk 



SQg 

disk 
l^dUs 



1 

-y^dUs 

./O 



{x' + A)- 



— cos(fcax' — 0o) 



^gUg 



ttA 27rA2 



S{\ko\aA) 



where 



S{\ko\aA) = IfcoctAI 



iw.{y)dy 



dx (37a) 
(37b) 

(38) 



,0 (y+ |A;o|aA)2- 

The first term in Eqn. fl37bp is the rate at which the undisturbed disk precesses the satel- 
lite's orbit. The second term is the rate at which the bending wave, whose amplitude is 
proportional to [ig by Eqn. (l20il . drives additional precession. The 5" function in that term 
is a dimensionless measure of the wave's contribution to the satellite's precession rate; that 
quantity depends on the wave's initial wavenumber |A;o|; and is plotted in Fig. [21 which shows 
that < S{\kQ\aA) < 1. 

Note that if the satellite instead orbited at the center of a narrow gap in the disk, 
then the first term in Eqn. fl37bp would be doubled due to the disk matter orbiting interior 
to the satellite. We might also expect additional precession to occur due to any bending 
waves launched in this interior disk, but it will be shown below that this contribution is 
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unimportant. With this in mind, we will generalize Eqn. fl37bl) to account for a possible 
inner disk by writing 



- - IS^(l^olaA) (39) 

disk TTA ZTtA^ 

where it is understood that e = 1 if the disk is entirely exterior to the satellite, and that 
e = 2 if the satellite instead orbits in the center of a gap whose fractional half-width is A. 
The satellite's total precession rate then becomes 



disk 



obi 



vrA 



27rA- 



.S(|A;o|aA) + tl. 



obi 



(40) 



When the disk and satellite are in steady-state, the satellite and its spiral bending pat- 
tern precess in concert, so Vtg = which after some manipulation yields the dispersion 
relation 

fJ'S 

Aid 

where 



nD{z)\k\aA = e + ^(l + ^) + 
jjLd V A/ 

/(|fcoaA|,^) = 5(|A;o|aA)- 



2/idA 
l^oflAp 



f{\kQaA\,z), 



(41) 



(42) 



is another function of distance z and wavenumber |/co|, one that is restricted to the interval 
-!</<!, and 



2l7r f RpA\ ^ 
yUc = — I — ^ 1 J2 



(43) 



which will be called the critical disk mass. 



2.2.3. limits on wave propagation 



The disk's ability to sustain these bending waves is assessed by multiplying the disper- 
sion relation (HTl) by v^f)/7rA, which yields 



HD 



V2I) 

ttA 



e+-{l + -r 
/id ^ A 



X 



+ 



t^sf 

2^idA 



(44) 



where the dimensionless wavenumber H = \^[)\k\a. Far downstream, where z ^ 1 and 



D(z) 



JDOO 



1 — e~ )/ H (see Eqns. [2^ and 155]) . so HD = 1 — e~ < 1. Wave propagation 
thus requires the right-hand side of Eqn. (jHl) to always be less than unity, which places an 
upper limit on the thickness of a disk that is able to sustain these bending waves, namely, 
that [) < [)max where 

7iA/V2 



£ + Hc/tld + l^s/'^tJ'd^'' 



(45) 
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upon setting a; = and / = 1 in order to obtain the most conservative limit on the disk's 
fractional thickness f)max- 

The remainder of this paper will assume that the disk is thin enough to sustain density 
waves, namely, that P) <^ f)max) or equivalently that if ^ 1, so that D ~ 0{1). Also recall 
that Section 12.21 anticipated a wave solution to occur when the satellite's mass is small. 
Specifically, when fis ^ SeyU^A, the right-most terms in Eqns. (14T1) or (l44l) may be neglected, 
which then provides the wavenumber A; as a simple function of distance x in the disk: 



where D{z) has been replaced with its average value over the first wavelength, D. And if the 
disk is sufficiently massive, namely, that fid ^ ^8fic/£, then the initial wavenumber at x = 
is |fco|ctA e/nD ~ 0.37, where e = 1 and D ~ 0.87 according to Fig. [2l In that limit, 
the first wavelength is Aq = 27r/|fco| — 2'K^DAa ~ 17A ■ a. However shorter wavelengths will 
result when /x^ does not exceeds the above threshold. 

Plugging Eqn. ( H6i) evaluated at x = into Eqn. ( l20i) then yields the wave amplitude in 
terms of the system's physical parameters: 



The satellite launches a spiral bending wave via its secular gravitational perturbations 
of the ring. Those perturbations tilt the orbital plane of the nearby ring particles, and they 
in turn tilt the orbits of the more distant parts of the disk. Tilting an annulus in the disk also 
tips that ring's angular momentum vector, so the excitation of a bending waves transmits 
in-plane angular momentum from the satellite to the disk. Consequently, wave-excitation 
damps the satellite's inclination Ig, and that rate can be calculated using the Lagrange 
planetary equation for Ig. 

The rate 6Is at which a single annulus in the disk damps the satellite's inclination Ig is 
Eqn. ([8]), again with n, a,I,Q —>■ Ug, ag, Ig, Qg and x' x' + A. Integrating the contributions 
by all annuli in the disk gives the satellite's total inclination-damping rate. 




(46) 



(47) 



Is 27rD/irfA 



2.2.4- inclination damping 




(48) 
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where ko is the initial wavenumber at the disk's inner edge, and the function 

cos{y)dy 



C{z) 



z + z 



cos 2 — Gliz) sinz 



(49) 



is shown in Fig. [31 Inserting Eqn. fl2Ul) into fH5]) and noting that /cq = — jA^ol then provides 
the inchnat ion-damping rate in terms of the system's physical parameters: 



C(|A;o|aA) /i^ 



2tt a 

The reciprocal of the above gives the e-fold timescale for the satellite's inclination decay: 



(50) 



orb 



C{\ko\aA)ij,s 



(51) 



where Porb = ^n/us is the satellite's orbit period. Note that for a planet that is not too oblate 
(e.g., Section EXID |A;o|aA = 0.37, so C(|A;o|aA) ~ 0.24 (Fig. [3]). This inclination-damping 
rate is also confirmed below, in a numerical simulation of spiral bending waves launched in 
a planetary ring. 



3. Simulations of spiral bending waves 



The rings model of lHahnI (120031 ) will be used to confirm the preceding results. The rings 
model treats the system as a set of N discrete gravitating annuli having semimajor axes aj, 
inclinations Ij, nodes Qj, and half-thicknesses hj. The model only considers the system's 
secular gravitational perturbations, so it also solves the same equations of motion, Eqns. 
([2]), but the model does so without making any of the wave-assumptions invoked in Section 
12.11 The model thus provides an independent check of the analytic results obtained above. 



3.1. waves in an exterior disk 

The rings model is used to simulate the spiral bending waves that are launched by 
an inclined satellite that orbits just interior to a disk. Figure S] shows the amplitude of 
this bending wave as it advances across a disk. The system's parameters are detailed in 
Fig. m Those parameters do not correspond to any real ring-satellite system; rather, these 
parameters were chosen to illustrate the results of Section [2] in the limit in which those 
results were obtained, namely, that the satellite's mass is small, i.e., fig <^ 2yUrfA so that 
Eqn. (H6|) is valid. Those parameters were also chosen so that the factor fic/ fJ'd appearing 
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in the wavenumber Eqn. fH6l) is 0.2, which causes the wavelength to slowly decrease with 
distance x as they propagate away. Nonetheless, the simulation reported in Fig. H] does 
correspond loosely to a small ~ 10 km satellite orbiting just interior to ring whose surface 
density is similar to Saturn's main A ring. 

Inspecting this system's angular momentum provides a quick check on the quality of 
this calculation. This system should conserve the in-plane component of its total angular 
moment um, Lj = ^ y^^mjUja'^Ij , where the sum runs over all rings and satellites in the 
system f lHahnll2003l ). The single-precision calculation shown in Fig. H] conserves L, with a 

-5 



fractional error of |ALj/Lj| < 2 x 10 

Note that the time for these waves to propagate a fractional radial distance x = Ar/a 

is 

_ Ar _ xPorb , . 

Cg 27r/irf 

where Cg is the waves' group velocity, Eqn. The simulated disk has a normahzed mass 

of /irf = 5 X 10~* and a fractional width x = 0.02, so the anticipated propagation time is 
tprop = 64 X 10^ orbits, which compares favorably with the simulation (see Fig. Hj). 

If the disturbance seen in Fig.lHis indeed a spiral bending wave, then the disk's longitude 
of ascending node Q{a) should steadily advance as a increases across the disk. This is 
confirmed in Fig. [5l which shows the waves' longitudes relative to the satellite's, Q{a) — Qg. 
These are also the longitudes where the warped disk passes through the central planet's 
equatorial plane. Note that this disk will have its maximum elevation at longitudes 90° 
ahead of that seen in Fig. [5l with its minimum elevation at longitudes 90° behind. Note 
also that the longitude of the disk's inner edge is 90° ahead of the satellite's longitude, as 
expected. And since the wavenumber k = —dfl/da is negative, this spiral pattern is indeed a 
leading wave. We also note that once the bending wave is established in the disk, the disk's 
longitudes precess at the same rate as the satellite's, i.e., Q{a) = Qg, and that the disk's 
inclinations are constant, /(a) = 0, which justifies our steady-state assumptions, Eqn. ([3]). 

Figure [5] also plots the dimensionless wavenumber |A;|aA across the disk at time t = 
75 X 10^. This is the moment when the bending wave is just starting to reflect at the disk's 
outer edge, which accounts for the curve's raggedness there. Also plotted is the expected 
wavenumber, Eqn. (146|) . which compares favorably. 

The rate at which the disk damps the simulated satellite's inclination Ig is shown in Fig. 
El where it is compared to the expected rate, Eqn. (15U]) . That rate is calculated by noting 
that the waves' initial wavenumber is I/cq^A ^ 0.63 at the disk's inner edge (see Fig. so 
the C that appears in Eqn. fl50l) is C{\ko\aA) = 0.32, according to Fig. [31 The expected and 
observed inclination damping rates are in good agreement. 
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3.2. Satellite in a gap 

The simulation described by Figs. H] and [5] is a bit of fiction, since there are no known 
satellites orbiting just interior to a broad planetary ring. For instance, all of the major 
Saturnian satellites orbit exterior to Saturn's main rings. However there are two noteworthy 
exceptions: the small satellite Pan, which orbits in the En cke gap in S aturn's A ring, and 



Daphnis, which inhabits the Keeler gap in Saturn's A ring (jPorc d 120051). 



rap 



A simulation of an inclined Pan as it orbits in the Encke gap is reported in Figure [TJ 
which shows the state of this system at time t = 7.5 x 10^ orbits. This is the time required 
for Pan to launch a leading spiral bending wave at the gap's outer edge that then propagates 
to the simulated ring-system's outer edge. Figure [7] shows that waves' initial wavelength 
is Aq = O.OOSTos ~ 500 km, and that their wavelength shrinks with distance x as they 
propagate towards the outer edge of the A ring, which lies a fractional distance x = 0.024 
away. So if these waves are not otherwise damped en route by collisions among ring particles, 
their wavenumber will have grown to \k\a ~ 1.1 x 10"^ when they reach the A ring's outer 
edge (see Eqn. H6|l . which corresponds to wavelength of A = 27r/|/c| ~ 80 km. 

Figure [8] also shows the system at the later time t = 1.5 x 10^ orbits, which is when the 
wave has since reflect at the simulated ring's outer edge and returned to the launch site. Here 
we see the superposition of an outbound leading wave with an inbound trailing wave, which 
results in a standing bending wave throughout the disk. As the Figure shows, when the 
standing wave emerges, it arranges the disk's longitudes fl such that they alternate between 
—90° and +90° of the satellite's node Qg at every half-wavelength. So if the waves launched 
by an inclined Pan do not get damped downstream, that bending wave will reflect at the 
outer A ring and return to Pan's vicinity where it can communicate its in-plane angular 
momentum back to the satellite. At this moment, inclination-damping then ceases. 

Figures [7] and [8] also show that Pan does not launch any inward-propagating waves. 
Although the satellite's secular perturbations do excite inclinations at the gap's inner edge, 
those disturbances do not travel further inwards. The wavenumber \k\ for any disturbance 
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that might propagate in the interior disk i: 



1 



) 



(53) 



k\ ~ 



TcDalA 



which is identical to Eqn. P6l) except for the sign on the disk mass term. Since the right 
hand side must be positive, Eqn. fl53|) tells us that waves in the inner disk can only propagate 
in the zone where < Xm, where 



is the distance of the waves' maximum excursion inwards of the satellite's orbit. Getting 
waves to propagate inwards a significant distance thus requires the disk mass to be sufficiently 
high, namely, fi^ ^ yUc/2, where e = 2 for a gap-embedded satellite. Saturn's A ring has a 
disk mass of /i^ ~ 5 x 10~^ and a critical disk mass of Hc = 7.8 x 10~^ (from Fig. [7] caption), 
so /id ^ /ic/2 is not well-satisfied, and inward-propagating bending waves are precluded. 

Also note that this simulation does not satisfy fid ^ ^8fic/£, which means that the 
wavenumber k does varies substantially across that first wavelength (see Eqn. H6l) . A 
wavenumber that is nearly constant over that first wavelength is of course a key assump- 
tion of Sections [2H3], so the analytic results obtained there might seem not apply to Pan. 
Nonetheless, when those formulas are compared to the model results, we find that Eqn. fH6l) 
to be in excellent agreement with the wavenumber k{x) exhibited by the simulated wave. 
But Eqn. fH7|) does overestimate the amplitude of this simulated wave by a factor of ~ 4, 
which in turn causes Eqn. ( 150|) to overestimate Pan's inclination-damping rate by the same 
factor. 

However Daphnis inhabits the narrower Keeler gap, whose fractional half width A = 
1.1 X 10~^ is a tenth that of the nearby Encke gap, so its /ic is 100 times smaller, and 
Xin — 120|A| = 0.013, which corresponds to a physical distance of x^ag ~ 1800 km, or 
about 10 wavelengths. So when the rings model is used to simulate the spiral waves that an 
inclined Daphnis would launch, we do indeed see a wave launched from the gap's inner edge. 



^This wavenumber is obtained by repeating the derivation of Section [2] by applying that to an annulus 
in the inner disk whose semimajor axis is again a = as(l + x + A), but with the distances x and A now 
understood as having negative values. So when integrating the net perturbation that the entire disk exerts 
on an annulus, the integration variable x' in Eqns. ^ and (|22p now runs over —oo < x' < ~x. The net effect 
of this is to merely change the sign on certain terms: Idisk is — IxEqn (jl3p . and Sk = sgn(fc) = sin^i = +1, 
where (f>i = 7r/2 is the longitude offset between the inner gap edge and the satellite. Any inward-propagating 
waves in the inner disk are trailing, since k > 0. The argument of the sinusoid in Eqn. (pSj) also changes 
sign. Accounting for these sign changes then yields Eqn. (|53p. 




(54) 
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That wave propagates inwards approximately a distance Xm, where it reflects and propagates 
outwards and across the Keeler gap. That satellite's simulated /-damping timescale is also 
in good agreement with the prediction, Eqn. fl5Tl) . 



4. External Vertical Resonances 

A satellite orbiting near a planetary ring also excites inclinations at its many external 
vertical resonances in the ring. This also communicates in-plane angular momentum be- 
tween the satell i te an d the ring, but in a manner that excites the satellite's inclination Ig. 
Borderies et al.l (Il984j ) calculate the rate at which the external resonances in a narrow ring 
of mass Sm' excite the satellite's inclination: 



gUsUslx 



M ' 



(55) 



where g = 0.0118, and x' is the satellite's fractional distance from the ring of mass Sm' = 
2Ti(Ta^dx' . The satellite's total /-excitation rate is the above integrated across the entire 
disk. Is = J^-^f^ 6 Is- If the satellite orbits in the center of a gap in a broad planetary ring, 
the total excitation rate due to external vertical resonances is 



A4 



(56) 



( iWard fc Hahrul2003l ). So if the satellite's orbit is to remain confined to the ring plane, this 
/-excitation due to the satellite's external resonances must be smaller than the /-damping 
that results from its secular interaction with the ring, Eqn. (l50l) . Comparing these two rates 
shows that the satellite's inclination is stable, i.e., Is < 0, when its gap is sufficiently large, 
namely, when 

2 ^ 27r5(/irf 



A' > 



(57) 



C(|A;o|aA)- 

Bending waves launched by Pan and Daphnis have initial wavenumbers of l/colctA 1 (see 
Eqn. H6l) . so C(|A;o|aA) ~ 0.3 (see Fig. [3]), and the requirement for inclination damping 
becomes A > 0.5y^. These satellites inhabit Saturn's A ring, which has a normalized disk 
mass of /irf ^ 5 X 10^*^, so their inclinations are stable if their gap half-widths are wider 
than A > 1.1 X 10~^. Pan easily satisfies this requirement (A = 0.0012), while Daphnis 
marginally so (A = 1.1 x 10~^). The /-damping timescale for Pan is quite short, only 
Ti ~ 1.7 X 10^ orbits (e.g., four times Eqn. ISTl) . which corresponds to ~ 2700 years. A 
comparable /-damp ing timescale is also obtained for Daphnis, whose size is about four times 
smaller than Pan's (ISpitale et al.ll2006l ). and whose gap is ten times narrower. 
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5. Summary and Conclusions 

The secular perturbations exerted by an inclined satellite orbiting in a gap in a broad 
planetary ring tends to excite the inclinations of the nearby ring particles. The ring's self 
gravity then allows that disturbance to radiate away in the form of a spiral bending wave. 
The wavelength A = 271 /\k\ of any outbound waves is obtained from Eqn. fH^ . which shows 
that A decreases with distance x from the nearby gap edge. These wavelength variations 
are due to a competition between the disk's self gravity and the differential precession that 
is due to the oblate central planet. As an example, we find that an inclined Pan, which 
inhabits the Encke gap in Saturn's main A ring, would excite a bending wave having an 
initial wavelength of about 500 km. And if that wave manages to propagate out to the outer 
edge of the A ring without damping, the wavelength will then have shrunk down to about 
80 km there. 

A gap-embedded satellite will also try to launch a wave at the gap's inner edge, but the 
range of these waves is limited by how far they can propagate until their wavenumber has 
shrunk to zero; see Eqn. (153|) . That distance is controlled by the width of the gap, with a 
narrower gap resulting in a greater inward excursion; see Eqns. (l54l) and (H3!) . For instance. 
Pan is unable to launch an inward-propagating wave, while Daphnis, which inhabits the 
narrower Keeler gap, could launch a disturbance that propagates inwards about 1% of its 
orbit before reflecting and propagating back out and across the Keeler gap. 

The amplitude of this wave is given by Eqn. (jlTj), and the excitation of this wave also 
damps the satellite's inclination quite vigorously, at a rate given by Eqn. fl50|) . This I 
damping mechanism also competes with the satellite's many other vertical re sonances in the 



ring, which try to pump up the satellite's inclination (iBorderies et al.lll984j ). However the 
secular /-damping will dominate when the satellite's gap is sufficiently wide, namely, when 
Eqn. (1571) is satisfied. Saturn's gap-embedded moon Pan satisfies this requirement, while 
Daphnis is at the threshold. This secular phenomenon also damps inclinations on a very short 
timescale r^, which for these satellites is of order 3000 yrs (Eqn. EI]). This of course assumes 
that these waves damp somewhere downstream, rather than refiecting at the ring's outer 
edge and returning to the launch site. But if these waves refiect and return without suffering 
significant damping, then a standing wave will emerge in the disk. That standing wave would 
also communicate some of its in-plane angular momentum back into the satellite's orbit, so 
further inclination-damping would cease. But if these waves instead damp downstream, 
then this secular phenomenon represents an important stabilizing infiuence that tends to 
confine a satellite's orbit to the ring plane. But this inclination-damping also shuts off any 
subsequent wave generation, so is seems unlikely that these waves might ever be observed in 
a planetary ring. 
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Finally, we note that the rings model employed here played a important role in guiding 
the analytic results obtained above. The model itself is a fairly easy-to-use set of IDL 
scripts, and other applications of this code are possible. For instance, the rings model has 
revealed that an eccentric satellite's secular perturbations can launch spiral density waves in 
a nearby ring, and the eccent ricity damp ing that is associated with that phenomena will be 



assessed in a foUowup study (IHahnl 120071 ). A copy of the rings model algorithm will also be 



made available to others upon request. 
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A. Appendix | 



The disk's vertical displacement is z{a,6,t) = asin/sin(^ — Q), where {a, 6) are the 
radial and azimuthal coordinates in the disk, and the inclination / and ascending node Q 
should be regarded as functions of distance a and time t. When a spiral bending wave is 
present in the disk, the longitudes have the form 



fl{a,t) = fi(ao,0) 



k{r)dr + ut, 



(Al) 



ao 



where fl{aQ, 0) is the longitude of the ascending node at some reference distance Oq at time 
t = 0, k{a) is the wavenumber of the spiral bending wave, and u is the angular rate at 
which the spiral pattern rotates, also known as the pattern speed. The signs in the above 
expression are chosen to follow the convention that a A; < spiral is a leading spiral, which 
means that a curve having a constant z{a, 6) in the disk traces a spiral that advances in 6 
as a increases. Thus the wavenumber k can be written as 



k = - 



and the group velocity is 



dQ 
da 

du 
dk 



(A2) 



(A3) 



(jToomrdll969l : iBinney fc Tremaind 119871 ) 
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Fig. 1. — A schematic showing the geometry of the ring-satelhte system, seen edge-on. A 
sateUite of mass and semimajor axis orbits interior to a broad planetary ring that 
extends to infinity. The sateUite's distance from the ring's inner edge is A in units of the 
sateUite's semimajor axis a^. A perturbed annulus in the ring has mass Sm, and it hes a 
fractional distance x away from the ring's inner edge, while the perturbing ring has mass 
Sm' and lies a fractional distance x' from the perturbed ring. 
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Fig. 2. — The functions z'^Ah{z) [from Eqn. IHM . solid curve], Bh{z) [Eqn. (Elj), dashed 
curve], D{z) [Eqn. (l33l) . sohd curve], and S{\kQ\a^) [Eqn. ( |38l) . dotted] are evaluated numer- 
ically for a thin disk having H = 0.01. These curves are plotted versus z = \k\ax, which is 
the dimensionless distance from the ring's inner edge in units of 27r wavelengths, or versus 
the dimensionless wavenumber |/co|cfA- The D{z) function is evaluated with |/c|aA = 0.37, 
a value that is justified in Section I2.2.2[ The average of D{z) over the first wavelength, 
< z < 27r, is -D ~ 0.87. Note also that z'^Ah{z) ~ cos 2; after the first wavelength. 
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Fig. 3. — The function C{\kQ\aA) from Eqn. fH9|) . plotted versus the dimensionless wavenum- 
ber |A;o|oA. 
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Fig. 4. — The rings model is used to simulate bending waves launched by an inclined 
satellite that orbits just interior to a disk. The satellite's normalized mass is /i^ = 10~^^, 
and the disk is comprised of = 500 rings having semimajor axes distributed over 1 + A < 
aj/tts < 1.02, where A = 5 x 10~^ is the fractional distance between the satellite and the 
innermost ring. The rings' fractional masses are /i^ = 3.9 x 10~^^, so the normalized disk 
mass is /x^ = nar'^/M = firf^lS/ag) = 5 x 10~^, where the rings' fractional separations are 
S/qs = 0.02/N = 4 X 10~^. The rings' fractional half-widths () is also set equal to their 
separations 6/as. The central planet's zonal harmonic is J2 = 0.012 and the planet's radius 
is Rp/tts = 0.45, so the system's critical disk mass is /ic = 1.0 x 10~* and l^'d = 0.2. 
The satellite's initial inclination is sinJ^ = 10~^, with all other rings initially having zero 
inclinations. The curves show the fractional amplitude of the bending wave, I{a)/Is, as it 
advances across the disk, shown at selected times t in units of 10^ orbital periods. The 
dashed hne is the expected wave amplitude, Eqn. fHTl) . with e = 1. 
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Fig. 5. — The upper figure shows the disk's longitude of ascending nodes Q{a) relative to the 
satellite's node fls, in units of vr, for the simulation of Fig. IHat time t = 75 x 10^ orbits, when 
the wave has swept across the disk. The lower figure shows the dimensionless wavenumber 
j/cjaA at this moment, where wavenumber is calculated from k = —dfl/da. Note that the 
simulated curve gets a bit ragged at the disk's outer edge, which is where the bending wave 
is just starting to reflect. The dashed line is the expected wavenumber, Eqn. (146!) . with 
e = 1. This spiral pattern has an initial wavenumber of about l^oloA = 0.63 at the disk's 
inner edge. 
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Fig. 6. — The rate at which the satelhte launching the wave in Fig. H] has its inchnation 
damped, Is, is plotted versus time t (in units of orbit periods). The solid gray curve is the 
expected rate, from Eqn. ( 150|) assuming |fco|oA = 0.63 and C{\ko\aA) = 0.32, where C is 
obtained from Fig. [31 
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Fig. 7. — A simulation of spiral bending waves launched by an inclined Pan orbiting in 
the Encke gap in Saturn's A ring. Figures show the disk inclinations I{a)/Is and relative 
longitudes Q{a) — Qg plotted versus semimajor axis a at time t = 7.5 x 10^ orbits, which is 
the time it takes the spiral bending wave to propagate across the simulated 'disk'. Note that 
the simulated disk is actually quite narrow since it only exte nds over 0.9 9 < a/a^ < 1.02. 
The system parameters are: Pan's mass fig = 8.7 X 10-12 Jporcol boosh . semimajor axis 
g.s = 1.34 X 10^ krn , inclination sini^ = 10"^, with an A ring surface density a = 50 gm/cm^ 
(IRosen et al.l Il99ll ) and a no rmalized disk mas s of /z^ = vraa^/M^ = 5 x 10^^. The gap 
half-width is Aa^ = 160 km (IBurns et al.ll2005l ). so its fractional half-width is A = 0.0012. 

= 500 rings are used to to simulate the wave in the disk exterior to the satellite, while 
50 rings are used to resolve the disturbance in the interior disk. The rings' fractional half- 
widths [) are set equal to their separations. Saturn's second zonal harmonic is J2 = 0.0163 
and the planet's radius is Rp = 0.45as, so the system's critical disk mass, Eqn. fH3l) . is 
/ie = 7.8 X 10-«. 
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Fig. 8. — The state of the system described in Fig. [7] at the lalfer time t = 1.5 x 10 orbits, 

which is when the bending wave has reflected at the simulated disk's outer edge and returned 

to the launch site, thereby estabhshing a standing wave in the disk. The fractional error in 

this system's total angular momentum is |ALj/Lj| < 5 x 10^^. 



